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igs: en eje de cuadratura de corriente estatorica. 


Resumen: El presente trabajo tiene como objetivo la descripcion del Va: Componente en eje directo de tensión rotórica. 
sistema de control del convertidor del lado rotor asociado a los e " en eje en cuadratura de tensión rotórica. 
generadores de inducción doblemente alimentados y su interacción iq " en eje directo de corriente rotórica. 


con el sistema de control de la turbina a efectos de lograr 
implementar el desacople de potencia activa y reactiva instantánea 
que fluye por el estator y la simulación de funcionamiento de la 


ig: " en eje de cuadratura de corriente rotórica. 
Wa: " en eje directo de flujo rotórico. 


máquina con este controlador, en caso que se encuentre trabajando Yar: " en eje en cuadratura de flujo rotórico. 
sobre una red de alta potencia de cortocircuito sin presentar Wa: " en eje directo de flujo estatórico. 
anomalías. Para dicho propósito se asumirá un modelo ideal de Ys: " en eje de cuadratura de flujo estatórico. 
convertidor (eficiencia unitaria y sin generación de armónicos), y del w,. frecuencia angular de la red eléctrica. 


modelo dinámico de quinto orden de la máquina de inducción, 
eligiendo un esquema de control vectorial donde el eje directo del 
marco de referencia dq es colineal con la tensión estatórica, girando 


w,: frecuencia angular eléctrica del rotor. 
@ş: velocidad de giro en el eje de alta velocidad. 


a la frecuencia de la red eléctrica. A partir del marco de referencia ES a desarrollado Boo el generador de inducción. 
rotatorio elegido se determina la interrelación existente entre las Pr: potencia activa instantanea por el rotor. 
diferentes magnitudes eléctricas involucradas en el análisis a efectos Ps: potencia activa instantánea por el estator. 
de deducir las expresiones que rigen el intercambio de potencia pi: potencia activa instantánea total. 

activa y reactiva instantáneas entre el estator y la red en función de qs: potencia reactiva instantánea por el estator. 
las corrientes en eje directo y en cuadratura del rotor, así como el Cp: coeficiente de potencia. 

acoplamiento existente de la dinámica de las componentes de la G: relación de la caja multiplicadora. 

tensión de rotor a través de estas corrientes. Esta última situación J: inercia del conjunto turbina — generador. 
conllevará el diseño de un controlador que desacople esta dinámica p: número de polos. 

mediante la utilización de lazos de alimentación directo 3 

(feedforward), permitiendo generar una referencia de tensión p: ángulo de paso. , , 

adecuada para que el convertidor del rotor pueda efectuar el A velocidad específica de la turbina. 
desacople de potencias. Se analizará la repercusión que tendrá esta v: velocidad del viento. 

forma de control sobre el desempeño global de la máquina sin R: radio de la turbina. 


detallar la intervención del convertidor de estator. Finalmente se 
realizarán distintas simulaciones a efectos de verificar el 
funcionamiento del generador con el controlador implementado. 


Palabras clave: Control de par, control de potencia reactiva, 
aerogeneradores de velocidad variable, generadores de inducción 
doblemente alimentados. 


I. NOMENCLATURA 


Vas. Componente en eje directo de tensión estatórica. 
vgs: " en eje en cuadratura de tensión estatórica. 





E ae e 
igs: " en eje directo de corriente estatorica. Fig.1: Esquema básico del DFIG y magnitudes eléctricas relevantes en el 


análisis. 
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| control del convertidor del rotor de un aerogenerador de 

velocidad variable se encuentra íntimamente ligado al 
sistema de control de la turbina: en función de las consignas de 
velocidad angular y par requeridas por la turbina para operar 
bajo una determinada premisa de control, el convertidor debe 
generar la tensión a inyectar en bornes del rotor de la máquina, 
con una amplitud y frecuencia adecuada a efectos de lograr los 
referidos requerimientos. En estos casos, es necesario que el 
generador trabaje sobre un amplio rango de velocidades, y a su 
vez, que desarrolle un par electromagnético acorde a la 
velocidad angular, a efectos de maximizar la extracción de 
potencia en el mayor rango de velocidad de viento posible, o 
limitar la misma, en función de las condicionantes físicas de 
los diferentes dispositivos que constituyen el aerogenerador. 


Paralelamente al control de par, lo que se traduce 
directamente en el control de potencia activa instantánea total 
que fluye desde generador hacia la red, es posible controlar en 
forma desacoplada, la potencia reactiva instantánea 
desarrollada por el estator. 


De acuerdo a las características de la fuente primaria, 
resulta imprescindible la utilización de un modelo dinámico de 
la máquina de inducción que permita obtener el grado de 
detalle necesario para observar la evolución instantánea de las 
diferentes magnitudes eléctricas que participan en el 
intercambio de potencia, y desarrollar una estrategia de control 
que se adapte a estas condiciones. Para abordar el análisis se 
ha optado por trabajar con el modelo de quinto orden de la 
máquina de inducción referida a ejes dq rotando a la 
frecuencia de la red a la cual se encuentra el generador 
interconectado y con el vector espacial de la tensión estatórica 
colineal con el eje d, de acuerdo a la figura 2. 


A Vs =v, + jO 


V 


S 


0(t) = o.t 


S 


» A 








Fig.2: Orientación de ejes dq elegidas para el análisis. 


Dado que el propósito principal de este trabajo es 
caracterizar el comportamiento dinámico del generador de 
inducción doblemente alimentado, la red elegida será de una 
alta potencia de cortocircuito y se supondrá que no habrán 
anomalías ni en la amplitud de tensión ni en la frecuencia de la 
misma. 


III. MODELO DINÁMICO DE LA MÁQUINA DE INDUCCIÓN 


De acuerdo a lo expuesto en la introducción, el modelo de 
la máquina de inducción utilizado será dado por las ecuaciones 


(1), donde explícitamente se está utilizando la orientación 
vectorial elegida (v,; = 0), y el par motriz de la turbina eólica: 


AY 

v, =ri, + S_@ l 
ds s“ds dt Was (1) 

AY, 
0 =r is += + OW a 
Var = rig ae (a, E o, w,, 

d 

v, =ri +r (0-0. 


qr r-qr dt 


3 
+ en Vai, g Y slas ) 


Y, = Lir + Ll y, 
Y, = Liy + bly 
War = Lig, t Lts 
VW. = Liy ELl 


El modelo presentado supone la máquina lineal, sin tomar 
en cuenta fenómenos como la saturación en el núcleo, efecto 
skin o armónicos espaciales. Por otra parte, la notación 
utilizada es de acuerdo a la convención de la máquina como 
motor. 


IV. IMPLEMENTACIÓN DEL CONTROL VECTORIAL 


A. Hipótesis de trabajo 


Si bien el modelo implementado por el simulador responde 
al modelo de quinto orden cuya descripción está dado por las 
ecuaciones (1), a efectos de realizar el control vectorial se 
realizarán las siguientes hipótesis de trabajo: 


1) La fuente a la cual se interconecta el aerogenerador 
es equilibrada, con amplitud y frecuencia constante. 
2) Las pérdidas por efecto Joule en el estator son 
despreciables: r, ~ 0. 
3) Los transitorios de flujo magnético en el estator son 
despreciables: Va = Vs =0 
dt dt 


La primera hipotesis de trabajo corresponde a la forma de 
operar el simulador; la segunda y tercer hipótesis serán 
verificadas mediante las simulaciones luego de efectuado el 
análisis que será desarrollado a continuación. En particular, la 
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tercera hipótesis presupone que el control se efectúa sobre una 
máquina de inducción de tercer orden. 


B. Condiciones particulares de trabajo 


Imponiendo las hipótesis planteadas en el conjunto de 
ecuaciones (1), se realiza una significativa simplificación en la 
ley de Faraday aplicada al estator del generador y a la segunda 
cardinal, que se terminan traduciendo en las siguientes 
expresiones: 


2 do, 1pARC,(4, Bj” 3p 
ne ne Pig O 
P —6G0, 
P 
Wis Lal + L nlar 7 0 
Wis = Ll a L mlgr m zi 


Se puede observar una simplificación importante en el par 
electromagnético desarrollado por el generador debido a la 
anulación de uno de sus términos. Por otra parte, se deduce 
que es posible controlar las corrientes en el estator si se logra 
controlar las corrientes que fluyen por el rotor de la máquina, 
lo cual se traducirá en el control de potencia activa y reactiva 
instantánea que fluye por el estator. Reordenando los términos, 
y recordando las expresiones de potencia activa y reactiva 
instantánea se deducen las siguientes ecuaciones: 


3p Goz 
nm Ve (3) 
Ps = aCe IF V oslas )= de ag 

2 De As 

3 VIV 
d, a Pato Vasa )= E a mig 


En síntesis, controlando la corriente en eje directo del rotor, 
se controla el par desarrollado por el aerogenerador, y 
controlando la corriente en eje de cuadratura del rotor se 
controlan los intercambios de potencia reactiva del estator. 
Asimismo, el control de estas corrientes implican la 
imposición de una tensión en el rotor del generador de 
inducción cuya amplitud y frecuencia se encuentra 
determinada por la dinámica de las dos componentes de la ley 
de Faraday del conjunto de ecuaciones (1). La particularidad 
de estas ecuaciones radica en el acoplamiento de las dinámicas 
de corriente y tensión del rotor en los respectivos ejes dq. Si 
no se actúa con el propósito de desacoplar la dinámica de las 
corrientes rotóricas, no será posible controlar en forma 
independiente el par y la potencia reactiva estatórica 
desarrollada por el generador de inducción. 


C. Desacople de dinámica de magnitudes del rotor 


Considerando las simplificaciones realizadas según el 
conjunto de ecuaciones (2), la ley de Faraday en el rotor de la 
máquina se reduce a las siguientes expresiones: 











di L 
a d ; 
Va Ag Oo q 0, 0,) di Vas 
ss 
di 
ene qr : 
Vor > Ploy +O dt +(0, — 0, Jai, (3) 


2 





L 
o=L.-—) 
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El acoplamiento referido en el conjunto de ecuaciones (3) 
implica que el controlador a diseñar debe considerar la 
inclusión de lazos de alimentación directa para eliminar la 
influencia de los términos cruzados, y asegurar un 
desempeño aceptable del control de cada una de las 
componentes de corriente rotóricas en lazo cerrado, de 
acuerdo al esquema de control que se observa en la Fig.3. 


(0,-0,). L_/(L@,) 





Fig.3: Controlador para el convertidor lado rotor de un DFIG 


Las referencias de las componentes de tensión rotórica en 
ejes dq generadas de esta forma permitirán actuar sobre el 
control PWM de forma tal que se implemente el control 
independiente de las corrientes en eje directo y en eje de 
cuadratura, por lo tanto, del par desarrollado, y la potencia 
reactiva intercambiada entre el estator de la máquina y la red 
eléctrica. Las simulaciones que serán presentadas presuponen 
que el ancho de banda del control PWM es muy alto y por lo 
tanto, la transferencia entre la tensión del convertidor rotórico 
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y la referencia generada es unitaria. 


V. ZONAS OPERATIVAS 


Si bien el objetivo principal del presente trabajo no es una 
descripción detallada del funcionamiento de las zonas 
operativas del generador, brevemente se hará una mención a 
las mismas debido a que son consecuencia directa de la 
implementación del sistema de control descripto 
anteriormente. 


Al requerir el control tanto de las componentes de corriente 
como de tensión en el rotor del generador de inducción, se 
controla la potencia que circulará por el rotor de la máquina, y 
al mismo tiempo, la potencia que circulará por el estator, 
motivo por el cual, el esquema de control presentado es 
relevante al punto que determinará la característica de los 
intercambios de potencia que tomarán lugar en el 
aerogenerador. 


Suponiendo que el aerogenerador se encuentra trabajando a 
una velocidad de viento constante, el modelo en régimen 
permanente derivará en las siguientes ecuaciones de 
intercambio de potencia (4), gráficamente representadas en la 
Fig.4: 


(4) 


Pa 2 3 PP. 
P, = 2 V gr l gy Y qr} (0, a o, o 


2 
m e s 


rap 3b 
P, Ps P, DE 








P,==Sp, 


Fig.4: Intercambios de potencia activa instantánea en un DFIG. 


Se puede observar que: 


1) En régimen de operación subsíncrona (s>0), el estator 
de la máquina generará potencia suficiente para 
abastecer las necesidades propias de consumo, 
determinadas por el sistema de control del 


convertidor estatórico, y entregar la potencia restante 
a la red eléctrica. 


2) En caso de operación sincrónica (s=0), no existiría, 
en términos ideales, intercambio de potencia activa a 
través del rotor de la máquina. 


3) En caso de operación supersincrónica (s>0), el flujo 
de potencia en el rotor se revierte, aportando potencia 
tanto rotor como estator a la red eléctrica. 


VI. SIMULACIONES 


Se presentarán dos escenarios de simulación en los cuales 
se mostrará la evolución de las corrientes desarrolladas en el 
rotor del generador de inducción frente a variaciones en la 
consigna de par y potencia reactiva en el estator. Al mismo 
tiempo se verificará la validez de las hipótesis realizadas 
para arribar a las ecuaciones que vinculan par y potencia 
reactiva con las respectivas componentes de corriente 
rotórica en eje directo y cuadratura. 


A. Simulación 1: Control de potencia reactiva instantánea a 


velocidad de viento constante. 


La simulación realizada consiste en mantener la 
velocidad de viento constante en un valor tal que la potencia 
total generada por la máquina es 750kW y variar las 
consignas de potencia reactiva desarrollada en el estator de 
la siguiente forma: inicialmente, el estator no intercambia 
potencia reactiva con la red; a los 2 segundos de simulación, 
se cambia la consigna de reactiva, pasando a producir 
150kVAr; a los 6 segundos, cambia nuevamente la consigna 
pasando a consumir SOKVAr, y finalmente a los 8 segundos 
se vuelve al estado inicial, es decir, no se intercambia 
potencia reactiva con la red (ver Fig.5). 


A partir de la simulación se puede observar que la 
evolución de las corrientes en eje de cuadratura de rotor y 
estator son las que regulan la evolución de la potencia 
reactiva, en tanto que las corrientes en eje directo, 
correspondientes al control de par (o potencia activa), 
permanecen aparentemente constantes (ver Fig.6 y 7) como 
era de esperar del análisis previamente realizado. Al realizar 
una ampliación de las corrientes en eje directo (ver Fig.8 y 
9), se puede observar la existencia de oscilaciones 
correspondientes a los transitorios de flujo magnético 
estatórico que fueron despreciados a efectos de llegar a las 
expresiones que vinculaban directamente al par y a la 
potencia reactiva a las corriente en eje directo y en 
cuadratura del rotor respectivamente, en los momentos que 
existe un cambio de consigna en la potencia reactiva a 
desarrollar por el estator. Esto indica que a pesar que el 
desacople no es perfecto, es razonable dada la amplitud y 
tiempo de asentamiento de los transitorios observados, 
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siendo estos mas rápidos que los transitorios de 
inicializacion del modelo dinamico (Fig.8 - Fig.9). 


Potencia activa y reactiva instantaneas (MW - MYAR) 

















Fig.5: Consignas de potencia activa y reactiva instantáneas. 
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Fig.6: Corrientes igr, Igs. 
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Fig.7: Corrientes las, lar. 


Corrientes en eje directo (4) 
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Fig.8: Ampliación ias 
Corrientes en eje directo (A) 


1062.2 


1062.1 




















Fig.9: Ampliación iar 


B. Simulación 2: Control de potencia reactiva instantánea 
frente a una rampa de velocidad de viento 


En esta simulación se conservará el patrón de consignas 
de potencia reactiva descriptas en la simulación anterior. La 
diferencia radica en que la velocidad del viento (según 
Fig.10) aumentará en forma de rampa pasando por las tres 
regiones operativas que definen el funcionamiento de la 
turbina eólica, estableciendo una consigna de par dinámica 
y variando la referencia de corriente en eje directa en el 
rotor, y por tanto en el estator (ver Fig.11). Inicialmente el 
generador entrega 750kW de potencia a la red y a los 0.5 
segundos de simulación comienza a aumentar la velocidad 
del viento, la evolución temporal de la potencia en estos 
instantes se debe a la aparición de la rampa en la zona de 
optimización de energía. A los 2.4 segundos entra en la 
zona de velocidades intermedias y a los 3.0 segundos se 
activa el control de ángulo de paso, limitando la potencia 
generada a 2.000kW. 
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Se puede observar en las Figs. 12 y 13 que la corriente en 
eje directo del rotor (y estator) controlan la dinamica de la 
potencia activa de la maquina, en tanto que las corrientes en 
eje de cuadratura controlan la dinámica de la potencia 
reactiva, de forma idéntica a lo expuesto en la primera 
simulación, observándose que es posible controlar en forma 
simultánea e independiente la potencia activa y la potencia 
reactiva. Finalmente en las Fig. 14, 15, 16, 17 se muestra el 
efecto del cambio de consigna de potencia reactiva y activa 
en las componentes de eje directo y de cuadratura de la 
corriente del rotor en forma respectiva, verificándose que 
los transitorios debidos a la componente transitoria del flujo 
del estator pueden ser despreciados y validando el esquema 
de desacople de potencia activa y reactiva propuesto. 


Velocidad del viento (m/s) 
15 T T T T T l 




















Fig.10: Velocidad de viento impuesta al aerogenerador 
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Fig.11: Consignas de potencia activa y reactiva instantáneas 
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Fig.12: Corrientes las, lar. 
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Fig.13: Corrientes igs, igr. 
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Fig.14: Detalle corriente ias frente a cambio de consigna de potencia 


reactiva. 
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Fig.15: Detalle corriente iar frente a cambio de consigna de potencia 
reactiva. 
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Fig.16: Detalle corriente ig, frente a cambio de zonas de operación de la 
turbina. 
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Fig.17: Detalle corriente ig, frente a cambio de zonas de operación de la 
turbina. 


VII. CONCLUSIONES 


A partir de la implementación del modelo detallado en el 
presente trabajo se concluye que es posible implementar el 
control desacoplado de potencia activa y reactiva en el 
generador de inducción doblemente alimentado tal como se 
realiza en las máquinas sincrónicas, a partir de la generación 
adecuada de las referencias de tensión y frecuencia para 
alimentar el rotor del generador, siendo importante observar 
que a diferencia del generador sincrónico, la 
implementación del control desacoplado de potencia activa 
y reactiva se implementa necesariamente con la máquina de 
inducción trabajando a velocidad variable. Las consignas de 
tensión y frecuencia en el rotor se imponen en función del 
par y velocidad de referencia que genera el sistema de 
control de la turbina, en tanto que la referencia de potencia 
reactiva que intercambia el estator y la red deberá ser fijada 
de acuerdo a la política de operación de la máquina. 


Por otra parte se debe recalcar la importancia de la 
validez de las hipótesis que conllevan a la reducción de 
orden del modelo de la máquina de inducción, a efectos de 
implementar el control, dado que permite visualizar 
claramente como se ven afectadas las diferentes magnitudes 
eléctricas de interés al actuar sobre las diferentes variables 
de estado. 


VIII. APÉNDICE 


A. Datos del generador 


U,: 690V — 50Hz 

Pa: 2.000kW 

p: 4 

Jy: 70 kgm? 

Ng max: 1620 r.p.m. 

rs: 1.5m0 

r’ 2.0mQ (referida al estator) 
Lint 2.4mH 

L,: 2.5mH 

L’ 2.5mH (referida al estator) 


B. Datos de la turbina y tren de engranajes 


La curva Cp(A,B) utilizada para las simulaciones 


correspondiente a la turbina de 2MW es la siguiente [2]: 





125 
C»(4,8)=0.22 "04-5 e ^ (7) 
l l 0.035 


A, A+0.086 B +1 

El maximo absoluto es Cp(6.325,0) = 0.4382 
P,: 2.000kW 
R=37.5m 
p = 1.225 kg/m 
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Ji: 1.4x10° kg.m? 
Mt max: 18 r.p.m. 
G = 90 
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